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ABSTRACT 

Context. The short burst GRB 130912A was detected by Swift, Fermi satellites, and several ground-based optical telescopes. Its X-ray 
light curve decayed with time normally. The optical emission, however, displayed a long-term plateau. 

Aims. We examine the physical origin of the X-ray and optical emission of short GRB 130912A. 

Methods. The afterglow emission was analysed and the light curve fitted numerically. 

Results. The canonical forward-shock model of the afterglow emission accounts for the X-ray and optical data self-consistently, so 
the energy injection model that has been widely adopted to interpret the shallowly decaying afterglow emission is not needed. 
Conclusions. The burst was born in a very-low density interstellar medium, which is consistent with the compact-object merger 
model. Significant amounts of the energy of the forward shock were given to accelerate the non-thermal electrons and amplify the 
magnetic fields (i.e., ~ 0.37 and eg ~ 0.16, respectively), which are much more than those inferred in most short-burst afterglow 

modelling and can explain why the long-lasting optical afterglow plateau is rare in short GRBs. 
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1. Introduction 


Gamma-ray bursts (GRBs) are short flashes of y-rays from 
the outer/deep space. Based on the duration distribution of the 
prompt emission, GRBs can be divided into two sub-groups 
dKouveliotou et al.1 fl993). One group has a typical duration of 
~ 20 s. The other has a much shorter duration < 2 s, cen¬ 
tered at ~ 0.1 s. The long GRBs are usually related to the 
death of massive stars, and one smoking-gun signature of such 
a scenario is the bright supernova emis_sion, as identified in 
most nearby long GRBs (IWooslevll993l:lMacFadven & Wooslevl 

1999: lHiorth et al.l2003l) . The short GRBs are rarer by a factor of 
about four according to the BATSE, or about ten from the Swift 
observation. The rate depends strongly on the energy bands and 
on sensitivities of the instruments (e.g., Qin et al. 2013, Zhang 
et al. 2012). The understanding of such a group of events had 
not been revolutionized until 2005 when Swift and HETE-ll lo¬ 
calized such events, and the long-l asting multi-wavelength after¬ 
glow emission had been detected (iGehrels et al.ll2005t iFox et aD 
2005). 


Although not as abundant as that of long GRBs, the after¬ 
glow emission data of short GRBs by then are valuable for re¬ 
vealing the physical processes taking place in the central engine. 
For example, the peculiar X-ray emission, such as the X-ray 
flares and the X-ray p l ateau followed by abrupt quick decline 
( Villasenor et alJ 120051: Barthelm v et al J 120051: IRowlinson et alJ 


1201 & Margutti et al. 20111) . was found to be inconsistent with 
the so-called standard forward-shock afterglow model. This 
emission ins tead implies the prolonged activity of the cen- 
tral engines dFan et alJl2005t iDai et alJfeoQd iPerna et alJl2006t 


Gao & Fanl 20061 Proga & Zhangl l2006t iMetzger et alJ 120081 : 
Rowlinson et al.l 20131) , which is possi bly associated with non- 


ignor able gravitational wave radiation dFan et al. 112013 a: Zhang 
12013 ). If these X-ray plateaus are indeed powered the inter¬ 
nal energy dissipation of supra-massive neutron star (SMNS) 
wind and the sharp declines mark the co llapse of the SMNSs 
dGao & Fanl 120061: IRowlinson etal.ll2013l) . then the maximum 
gravitational m ass of non-rotating neutron star c an be estimated 
to be ~ 2.3 M 0 dFan et al. ll2013atlLi et al. Il2014l) . Moreover, the 
detection of a weak infrared bump in GRB 130603B strongly fa- 
vors the physical or i gin of the merger o f a compact ob ject binary 
dTanvir et al J 120 1 3l Berger etaD 12013b I.Tin et alJl2013l) . To ac¬ 


count for the shallowly decaying X-ray emission, a millisecond 
magnetar central engine, which was still active in - 10 3 s after 


double neutron stars (Fa 

n et al. 2013b; de Ugarte Postigo et al. 

2014j; Lii & Zhang 20 5 

Lii et aLll2015l). 


Motivated by this progress, in this work we examine the 
physical origin of the X-ray and optical emission of short 
GRB 130912A, which is characterized by a long-lasting opti¬ 
cal plateau. In section 2 we introduce the observations of GRB 
130912A. In section 3 we interpret the data. We summarize our 
results with some discussion in section 4. 


2. Observations 

2.1. Swift observations 

The Swift Burst Alert Telescope (BAT) detected GRB 130912A 
at 08:34:57 UT on September 12, 2013 dD’Eliaet al.ll2013i) . The 
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7^90 duration is 0.28 + 0.03 second. The light curve shows two 
overlapping peaks. The time-averaged spectrum of the first 0.32 
second is best fitted by a simple power law, with photon index 
r r = 1.20 ± 0.20. The total fluence is 1.7 + 0.2 x 10~ 7 erg cm 
2 , and the peak photon flux is 2.2 + 0.3 ph cm" 2 s~ 1 (Krimm et 
al. 2013). All values are in the 15 - 150 keV energy band. There 
is no evidence of extended emission detected in the BAT energy 
range dNorris et al.l2010l) . which makes it an unambiguous short 
GRB. 

The Swift X-ray Telescope (XRT) began to observe the field 
at 08:36:31.7 UT, 93.9 seconds after the BAT trigger. Kennea 
et al. (2013) analyzed the initial XRT data and report that the 
light curve can be modelled with a power-law decay with a de¬ 
cay index of a = 1.20 + 0.04. The spectrum formed from the PC 
mode data can be fitted with an absorbed power law with a pho¬ 
ton spectral index of F x = 1 .57/[j‘ 2 g and Ah = 1.49^° ® x 10 21 
cm -2 , consistent with the Galactic value of Ah = 1 . 2 x 1 0 21 cm" 2 
(Kalberla et al. 2005). 

The Swift UVOT took 150 seconds to find the chart exposure 
98 seconds after the BAT trigger dD’Elia et al.ll201 3l). No optical 
afterglow within the enhanced XRT position (iBeardmore et alj 
1201 3li h as been detected in this initial and the subsequence ex¬ 
posures ((Chester & D’Ehall^Oll . 

2.2. Ground-based optical observations 

The field of GRB 130912A was observed at early times by sev¬ 
eral ground-based optical telescopes, but only two detected the 
afterglow. GROND started observing at 08:50 UT on September 
12, 2013 (16 min after the GRB trigger), and found the af¬ 
terglow at coordinates RA (J2000) = 03:10:22.23 and Dec. 
(J2000) = 13:59:48.7 within the Swift XRT error box. Over 
the first hour-long observation in the r' filter, the afterglow 
seemed to be constant, with a magnitude (here and through¬ 
out thi s pape r , ma gnitudes are in AB system) of r’ =22.0+0.2 
(iTanga et al.ll2013l) . Another telescope, P60, started observing at 
8:48 UT on September 12, 2013 (13 min after the Swift trig¬ 
ger) also in r' band, confirmed the GROND observation that 
the source was almost unchanged in the first hour. It led to the 
suspicion that the detected source is the host galaxy of GRB 
130912A ( Cenko e t al . 2013 1: however, when RATIR observed 
the field again started on 2013 September 13.30, the source had 
faded down to a limited magnitude 23.89 (3-sigma) in SDSS 
r band, so the previously detected source should be the after¬ 
glow (iButler et al.ll2013l) . According to the RATIR observation, 
the host galaxy is fainter than 23.89, 23.79, 22.86, 22.38, 22.30, 
and 2T78 magnitudes in r, i, Z, Y. J, and //bands, respectively 
(iButler et all2013l) . None of these ground-based observations is 
publicly available in the form of raw data. We collected the pub¬ 
lic GCN data as shown in Table 1, and loose limits had been 
discarded. The magnitudes used in the following analysis were 
corrected for the Galactic extinction, assuming E(B -V) = 0.28 
dSchlegel et al. 1 998 ) and a ratio of total-to-selective extinction 
Rv = 3.1, the Galactic extinction in r' band is A r > = 0.78. 


3. Interpreting the optical and X-ray afterglow 

In the Swifts ra, long-lasting plateau-like X-ray emission (i.e., 
the so-called shallow decay phase) was detected well in a good 
fraction of GRB afterglows (e.g., Zhang et al. 2006; Nousek et 
al. 2006). The leading interpretation of the long-lasting plateau¬ 
like X-ray emission is the energy injection model, which is valid 
if the central engine works continually, or alternatively, the bulk 


Lorentz factor of the outflow material has a wide distribution 
(e.g., Dai & Lu 1998; Zhang & Meszaros 2001; Fan & Xu 
2006). A general prediction of the energy injection model is that 
the temporal behaviours of multi-wavelength afterglow emission 
will be shaped simultaneously (Fan & Piran 2006). However, the 
X-ray and optical observations of the GRB afterglow usually do 
not track each other. The X-ray afterglows usually displayed an 
early shallow decline, which is often not observed in the opti¬ 
cal (e.g.. Fan & Piran 2006; Panaitescu et al. 2006). The phys¬ 
ical reason is still not clear for such a puzzle, for example, the 
widely-adopted energy injection model for the X-ray decline is 
found to be unable to account for the optical data. 

Recently, there is another unusual situation. As shown in 
Fig.l, the X-ray afterglow emission of GRB 130912A can be 
fitted by a single power law (Beardmore et al. 2013), as found 
in most GRB X-ray emission. However, the optical emission is 
plateau-like on a very long timescale ~ 3.2 x 10 3 s, which is a 
very strange behaviour that is rarely observed in the optical af¬ 
terglow. A long-lasting optical plateau was also detected in GRB 
090510 (De Pasquale et al. 2010; Gao et al. 2009). The duration 
of the optical plateau of GRB 090510, however, is just about half 
of that of GRB 130912A. The optical plateau of GRB 130912A 
is thus very likely the longest one people have detected so far in 
short GRB afterglows. In reality, the lack of optical variability 
in the first hour after the trigger of GRB 130912A motivated the 
idea that these emission were from the host galaxy (Tanga et al. 
2013; Cenko et al. 2013). The afterglow emission nature of the 
optical data had not been established until significant fading was 
identified about one day after the burst (Butler et al. 2013). The 
main purpose of this work is to interpret the “unusual” optical 
plateau and the X-ray emission self-consistently. 

Because nothing is unusual in the X-ray band for GRB 
130912A, we conclude that the energy injection model does not 
work for the current data, as discussed above. It is widely known 
that the forward-shock emission is governed by some physical 
parameters that can be parameterized as (e.g.. Sari et al. 1998; 
Yost et al. 2003, Fan & Piran 2006) 

F y , m ax = 6.6 mJy i^-^Dl^ M e l J^ 2 E K53 n 0 12 , (1) 


v m = 2.4 x 10 16 Hz E l ^ 2 _ 2 el_ 


v c = 4.4 x 10 16 Hz E k \ /2 e B 3 ! 2 n Q 


(l+ZY 112 -112 _1_ 

V 2 / rf -3(i + T) 2 ’ 


( 2 ) 

(3) 


where C p = 13 (p - 2)/[3 (p - 1)], (eg) is the fraction of shock 
energy given to the electrons (magnetic field), tj is the time in 
days, the Compton parameter Y ~ (-1 + 1 + 4Tje e leB)l2, rj ~ 

minjl, (v m /v c ) (p ~ 2 ^ 2 ), and v c = (1 + Y) 2 v c . Here and throughout 
this text, the convention Q x = Q/10 x has been adopted. 

At t ~ 3.2 x 10 3 s, if we have min{v c , v m } = v c > 5 x 10 14 Hz, 
the temporal behaviour of the optical emission would be F v oc 

FymaxVc oc t 1 ! 6 . and the X-ray emission light curve should be 
F Vx oc C^maxVc' 2 oc t~ in the case of v c < v x < v m or, al¬ 
ternatively, F Vx oc Fy maxVc^Vm '^ 2 00 t~ (2p ~ 2 ^ 4 in the case of 
Vc < Vm < Vx- While the temporal behaviour of optical emission 
agrees nicely with the data, the X-ray emission does not. The 
case of v c < v x < v m is clearly at odds with the data. The case 
of v c < v m < v x on the timescale of 100 s < t < 3.2 x 10 3 s 
is also inconsistent with the X-ray spectrum F v oc y -o.50±o.i6 
(http://www.swift.ac. uk/xrt _spec tra/). 
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Table 1. Optical observations of the field of GRB 130912A 


Telescope 

Data start 
(UT) 

Observation time after trigger 
(s) 

Filter 3 

Magnitude 11 

Flux 

(erg cm^s^'Hz -1 ) 

GROND 

08:50/12/09/2013 

960 

GROND r’ 

22.0 ± 0.20 U) 

5.75 x 10^ 

P60 

08:58/12/09/2013 

1380 

GROND r’ 

21.77 ± 0.20 (2) 

7.11 x 10~ 29 

P60 

09:28/12/09/2013 

3180 

GROND r’ 

22.09 ± 0.25 (2) 

5.29 x 10^ 29 

RATIR 

07:33/13/09/2013 

~ 8.1 x 10 4 

SDSS r 

> 23.89 {3) 

< 1.01 x 10- 29 


a. The difference between GROND r' and SDSS r is less than 0.04 mag assuming a power-law spectrum with an index /? between 1 and -2. Here 
the flux of the afterglow is expressed as F v oc fv 13 , 

b. The flux is reported with the ltr statistical error, and the upper limit is at the confidence level of 3<r. References: (l1 Tanga et al. (2013), (2) Cenko 
et al. (2013), (3) Butler et al. (2013). 



Fig. 1. X-ray (square) and optical (circle) light curves of GRB 130912A and the theoretical model curves. The X-ray data are from 
the UK Swift Science Data Center (Evans et al. 2009) and transformed to 1.732 keV. In order to reduce the influence of the error 
of the spectral index on the flux calculation, we used the geometric mean of the lower and upper boundaries of the corresponding 
X-ray energy band 0.3 - 10 keV. Proper corrections for extinction in the Milky Way Galaxy have been made. The solid and dashed 
curves are the theoretical optical and X-ray afterglow prediction with a forward shock. 


If we, instead, have min{v c ,v m } = v m > 5 X 10 14 Hz at 
t ~ 3.2 x 10* s, the temporal behaviour of the optical emis¬ 
sion would be F Vopl oc F v , max Vm^ 3 °c fO 2 , and the X-ray emission 
light curve should be F Vx oc F vmM v[y '^ 2 oc /-30-0/4 in the case 
of v m < v x < v c or, alternatively, F Vx oc F Vjmax v,J /2 v^ 1|/2 oc 
f -(3p-2)/4 j n t h e case 0 f y m < v c < v x . Now the temporal be¬ 
haviours of both optical and X-ray emission are consistent with 
the data for p ~ 2.3, as are the spectral behaviours. Since an op¬ 
tical plateau lasting a few thousand seconds is rare in the short 
GRB afterglow data, it is highly necessary to examine whether 
the required forward-shock physical parameters are reasonable 
or not. 


To self-consistently interpret the optical and X-ray data, 
v c (f ~ 1.0xl0 4 s) «;7.254xl0 lf) Hz, v m (f ~ 3.2xl0 3 s) « 5xl0 14 
Hz and F Vjmax ~ 0.02 mJy are needed. Following Zhang et al. 
(2015), we have 

e B.-2 Ek -^ n 0 2 ~ a ’ 

(4) 

<5 2 3^2<-l ^ / 2 - 

(5) 

+ Y)~ 2 - O 

(6) 


where 


—F n 

6.6 v " max L,28 


10- I6 v m c; 2 (y 


- 1/2 


M' 


b = — 

u 2.4 


f^ 2 3 , and c = 


x KT lf V ( 


\ 1/2 


1/2 


(l±z) 

4.4 - 2 J V.-3- 

Now we have three relations but four free parameters (i.e., 
Ek, f b, no), which implies that these parameters cannot be 
uniquely determined. However, (Ek, e e , ee) can be expressed as 
the functions of no, and it is possible to reasonably constrain the 
range of «o, as found in Zhang et al. (2015). 

In the case of Y < 1 (i.e., the synchrotron-self Compton cool¬ 
ing is unimportant), (1 + Y)~ 2 can be ignored, and we have (see 
also Zhang et al. 2015) 


€b- 2 = a 3 c 5 n 


(7) 


-i = a 5 fo 5 c io ? n 


Ek, 53 = <UcUi 0 


( 8 ) 

(9) 


These three variables weakly depend on hq. 

In the case of Y > 1 (i.e., the synchrotron-self Compton ra¬ 
diation is important), we have (1 + Y) 2 a T~ 2 and then 


_S_ 2[p-l) 4(p-l) . 

2 = (10cc/) 5 ?- 9 a 5 p- 9 b 5 p- 9 n. 


3/2+1 
" 5/7—9 


( 10 ) 
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1 P-2 3p-5 

€ e -1 = (10cJ) 5 ^ 9 fl 5/7-9 £ 5/7-9 ft 


P-1 

5/7-9 

0 


( 11 ) 


4 6y-10 2(p-l) -i-l 

^a -,53 = (10«/) 5 p- 9 a 5 p- 9 b 5 p- 9 n 0 F , 


( 12 ) 


P~± 

where d = 2 . Compared with the case of Y < 

1, 6b and Ek depend strongly on hq. 

The redshift z is unknown, and we assume that it is in the 
reasonable range of 0.1 < z < 1.4. However, for z < 0.4 in 
solving eqs. dm we do not find any reasonable values of the 
parameters with any given no, while solutions are obtainable for 
largerz. In Fig. [2]we present the physical parameters (£L e e , fo) 
as functions of no forz = (0.4, 0.6, 0.8, 1.0, 1.2, 1.4), respec¬ 
tively. As already shown in the analytical approaches, both E^ 
and 6 b evolve with no quickly, while the dependence of e L . on 
no is rather weak. One can also find from the figure that there 
is a very interesting constraint that no < 3 x 10~ 3 crrT 3 for a 
reasonable e e < 0.4 (i.e., expected to be not much larger than 
the equipartition value ~ 1/3), though the afterglow physical 
parameters cannot be uniquely determined. GRB 130912A was 
therefore born in a very low-density medium, consistent with the 
compact-object merger model. 

To better show whether the afterglow model can indeed rea¬ 
sonably account for the data or not, in Fig JT] we numerically 
fit the optical and X-ray data of GRB 130912A. The numer¬ 
ical calculation code was developed by Fan & Piran (2006) 
and IZhang et~afl (120061) . In it, (i) the dynamica l evolution of 
the outflow formulated by iHuang et al.1 (l2000h that can de¬ 
scribe the dynamical evolution of the outflow for both the rel¬ 
ativistic and non-relativistic phases has been adopted, (ii) The 
energy distribution of the shock-accelerated electrons is cal¬ 
culated by solving the continuity equation with the power- 
law source function Q oc y~ p , normalized by a local injec¬ 
tion rate (iModerski et al.l l2000t) . (iii) The cooling of the elec¬ 
trons due to both synchrotron and synchrotron-self Compton 
has been taken into account. Assuming z = 0.72 (following 
iRowlinson et al.l (1201 3l) we adopt the average value of redshift 
of short GRBs), the fit parameters are (E^, no, e c , e b, p, 0 t ) = 
(1.7x 10 51 erg, 0.002, 0.37, 0.16, 2.3, 0.03), where <9j is the half 
opening angle of the GRB ejecta. An isotropic fireball is found 
to be unable to reproduce the data. The inferred e e and eg are at 
the high end of the distribution of the shock parameters of short 
GRB afterglows (e.g., Soderberg et al. 2006; De Pasquale et al. 
2010), which is unexpected since the optical afterglow plateau 
of GRB 130912A has the longest duration people have ever de¬ 
tected in short events. We would like to point out that the above 
fit parameters are for A = 0.78, which is very high. If A r > is 
intrinsically smaller, F Vmax and hence a are lowered accordingly. 
As shown in eqs.(7-9), or alternatively eqs.(10-12), eB and e e 
would increase, while would decrease. The contrary holds in 
the case of larger A r >. 


4. Summary and conclusions 

The most remarkable feature of the short burst GRB 130912A is 
an optical plateau lasting about 4000 s, which is the longest one 
in current short GRB observations, and it is about twice longer 
than that of GRB 090510. In this work we examined whether any 
’’unusual” information can be extracted from the afterglow data 
of GRB 130912A. Though the energy injection model has been 
widely adopted to interpret the shallowly decaying afterglow 
emission of long and short GRBs (see Zhang et al. 2006; Nousek 
et al. 2006 and the references therein), it was found to be un¬ 
able to account for the X-ray and optical data of GRB 130912A 


self-consistently. Instead the canonical afterglow emission of an 
ejecta with an opening angle (9, ~ 0.03 can reasonably reproduce 
the data. The circum-burst medium is found to be ISM-like and 
has a very low density ~ 10~ 3 cm -3 , consistent with the model 
of merger of binary compact objects (either double neutron stars 
or a neutron star-black hole). Significant amounts of the energy 
of the forward shock were given to accelerate the non-thermal 
electrons and amplify the magnetic fields (i.e., e c ~ 0.37 and 
6 b ~ 0.16, respectively), which are much more than those in¬ 
ferred in most short burst afterglow modelling and which can 
explain why the long-lasting optical afterglow plateau is rare in 
short GRBs. 
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